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Abstract. Extremely Red Objects (EROs, R—K>5) constitute a heterogenous class of extragalactic sources including high 
redshift elliptical galaxies, dusty star-forming systems and heavily obscured AGNs. Hard X-ray observations provide an 
unique and powerful tool to uncover obscured nuclear or star-forming activity. We present the results of XMM-Newton 
observations of the largest sample of EROs available to date (about 450 objects over a contiguous area of 700 arcmin 2 ). 
Five of the 46 hard X-ray selected sources brighter than 3x 10~ 15 cgs in the 2-10 keV band, are associated with EROs. All 
of the X-ray detected EROs show rather extreme X-ray-to-optical flux ratios, suggesting the presence of highly obscured 
AGN activity. We also report on the X-ray stacking analysis of spectroscopically identified EROs in the Chandra Deep Field 
South. 
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1. Introduction 

Extremely Red Objects (EROs, R — K> 5) show the bulk of 
the emission in the near-infrared band and are associated with 
extremely faint optical counterparts that have so far limited 
the identification process. The first available optical-infrared 
spectra (e.g., Cimatti et al. 1998; Smail et al. 1999; Pierre 
et al. 2001) and the recent results from the K2 0 survey 
(Cimatti et al. 2002) indicate that EROs counterparts lie at 
high redshifts. In particular, their colors are consistent with 
those of early-type passively evolving galaxies, dusty star¬ 
forming systems and Active Galactic Nuclei (AGN) reddened 
by dust and gas in the redshift range z = 0.8-7 3. The relative 
fraction of these different subclasses among the ERO popu¬ 
lation is a key parameter in the study of the galaxy evolution 
and can constrain models which link the formation of massive 
elliptical galaxies and the onset of AGN activity (Granato et 
al. 2001; Almaini et al. 2002). 

Hard X-ray observations provide a powerful tool to uncover 
AGN among the EROs population. Indeed, a sizeable fraction 
of hard X-ray sources recently discovered in deep Chandra 
and XMM-Newton surveys are associated with EROs, the 
exact value depending on the limiting fluxes reached in the 
optical and X-ray bands (Brusa et al. 2002a; Mainieri et al. 
2002; Alexander et al. 2002). Moreover, the excellent imag¬ 
ing capabilities of Chandra enable to probe the average X- 
ray properties of EROs even beyond the limiting flux of deep 
surveys, using the stacking analysis technique. 
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In this framework we have developed two complementary 
programs of hard X-ray observations of EROs making use 
of both Chandra and XMM -Newton data. First, we present 
the results of a XMM -Newton survey on a large, complete 
sample of EROs which is the most suitable one to estimate 
the fraction of AGN-powered EROs among the optically se¬ 
lected population. Then, we have constrained the high-energy 
properties of a sample of spectroscopically identified non- 
AGN EROs in the Chandra Deep Field South (CDFS) via 
stacking analysis. 

Throughout the paper, a cosmology with Hq = 70 km s -1 
Mpc^ 1 , Cl m = 0-3 and Oa=0.7 is adopted. 

2. XMM-Newton observation of the “Daddi 
Field” 

2.1. X-ray data analysis and results 

We have started an extensive program of multiwavelength ob¬ 
servations of the largest sample of EROs available to date 
(~450 sources), selected in a contiguous area over a ~ 700 
arcmin 2 field (the “Daddi field”, Daddi et al. 2000) and com¬ 
plete to a magnitude limit of K=19.2. Deep optical (R-~ 26.2 
at the 3cr level) photometry is available (see Daddi et al. 2000 
for details in optical and near-infrared data reduction) and 
VIMOS spectroscopy is planned at VLT. 

The Daddi field was observed by XMM -Newton on August 

3, 2001 for a nominal exposure time of ~ 50 ks. The XMM- 
Newton data were processed using version 5.2 of the Sci¬ 
ence Analysis System (SAS). The event files were cleaned 
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Fig. 1. R-band magnitude vs. the 0.5-8 keV flux for the 91 sources 
detected in the XMM -Newton observation. Squares are sources de¬ 
tected both in the K-band (K< 19.2) and in the R-band (R< 26.2) 
images, triangles are sources detected in the R-band image only. En¬ 
larged squares are EROs. The dashed area corresponds to the region 
of 0.1 <iji a < 10. 


up from hot pixels and soft proton flares; the resulting ex¬ 
posure time is ~ 35 ks in the MOS1 and MOS2 detectors 
and ~ 30 ks in the pn detector. The EBOXDETECT task, 
the standard SAS sliding box cell detect algorithm, was run 
on the 0.5-2, 2-8 and 0.5-8 keV cleaned events. We limited 
the X-ray analysis in a region of ten arcmin radius from the 
centre of the XMM-Newton pointing, to take into account 
the broadening of the XMM -Newton PSF at increasing dis¬ 
tance from the aim point. Ninety-one sources were detected 
in the full X-ray band with a threshold in Poisson probability 
of 5xl0 -6 ; the corrisponding flux limit is ~ 3 x 10 -15 erg 
cm -2 s -1 assuming a power-law spectrum with I’= 1.8 and 
Galactic absorption (N h = 5 x 10 20 cm -2 ). 

The X-ray centroids have been corrected for systematic er¬ 
rors with respect to the optical positions of three bright 
quasars in the field (Hall et al. 2000) and then were 
cross-correlated with the /T-barul and /(-band catalogs. We 
searched for near-infrared counterparts within a radius of 6" 
from the X-ray position; we note, however, that the 75% of 
the sources lie within a 3" radius circle. The results and the 
optical identifications are summarized in Fig. 1: 80% of the 
0.5-8 keV sources have at least one near-infrared counterpart 
down to the K-band magnitude limit; four sources are not de¬ 
tected at the R-band limit. Nine out of ninety-one sources are 
associated with EROs (~ 10%); the fraction is the same (5/46) 
for the 2-8 keV sample. We refer to Brusa et al. (2002a) for a 
more detailed discussion on the XMM -Newton data analysis, 
the associations of X-ray sources with optical counterparts 
and the analysis of X-ray colors. 


2.2. AGN fraction in optically selected samples 

The AGN fraction among optically selected EROs samples 
is still very uncertain (Brusa et al. 2002a; Mainieri et al. 
2002), and strongly dependent on the limiting fluxes reached 
in both the X-ray and near-infrared bands. About 370 EROs 
fall within the XMM -Newton area analysed in this work. The 
fraction of AGN-powered, hard X-ray selected EROs at the 
near-infrared (K=19.2) and X-ray flux limits of our survey is 
therefore ~ lq-2% (5/370). This fraction rises up to 14% at at 
K<20.1 and S2 -io~ 4 x 10 -16 erg cm -2 s -1 limiting fluxes 
of the CDFN observation (Alexander et al. 2002), supporting 
the idea that the bulk of the EROs population is not related to 
active phenomena. 

It is worth noticing that these AGN-powered EROs, as well 
as others examples of hard X-ray selected EROs show a rela¬ 
tively well-defined correlation between the X-ray fluxes and 
the optical magnitude around fx/fo — 10. This correlation 
is shifted by one order of magnitude from the one found by 
ROSAT for soft X-ray selected quasars (Hasinger et al. 1998) 
and recently extended by Chandra and XMM-Newton obser¬ 
vations to hard X-ray selected sources (Alexander et al. 2002; 
Lehmann et al. 2001). The most plausible explanation of such 
a high X-ray-to-optical flux ratio in these objects is the pres¬ 
ence of high obscuration towards the active nuclear source. 
This hypothesis is strongly supported by optical identifica¬ 
tions available for a handful of objects (Cowie et al. 2001; 
Mainieri et al. 2002; Hornschemeier et al. 2001) and by the 
X-ray spectral analysis (see Sect. 4). 

3. Stacking EROs in the CDFS 

The K2 0 EROs sample (Cimatti et al. 2002) in the Chan¬ 
dra Deep Field South (CDFS) area includes 48 objects at 
K s < 20. For the purposes of the present work we con¬ 
sider only the 21 spectroscopically identified EROs which 
have been classified in two categories: 13 “dusty”, i.e. objects 
showing at least an emission line (tipically [Oil] A3727) over 
a red continuum, and 8 “old”, i.e. objects with an absorption¬ 
line spectrum, consistent with that of early-type, passively 
evolving galaxies, distributed on a similar range of redshifts 
(z = 0.8 -i-1.6). 

For the X-ray analysis, we made use of the public CDFS one 
megasecond Chandra observation available in the archive 
(Giacconi et al. 2002). At first, we have searched for in¬ 
dividual X-ray emission from all the objects in the EROs 
spectroscopic redshift catalog of Cimatti et al. (2002). Only 
one source, a dusty ERO at z = 1.327, was individually 
detected in the one megasecond exposure. The optical po¬ 
sition is almost coincident (A < 0.3”, within the accu¬ 
racy of the K-band position) with the Chandra source CXO 
CDFS J033213.9-274526 in the Giacconi et al. (2002) cata¬ 
log. The very hard X-ray spectrum and the high intrinsic X- 
ray luminosity unambiguously reveal the presence of an ob¬ 
scured AGN in this source. Indeed, this object is undetected 
in the soft band, implying an intrinsic column density larger 
than 4xl0 23 cm -2 , and an unabsorbed hard X-ray luminos¬ 
ity L 2 _io > 4 x 10 43 erg s -1 ; its properties are then more 
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Fig. 2. Stacked full-band images of the “dusty” (left panel ) and 
“old” (right panel) EROs. The images are 9x9 arcsec and have been 
smoothed with a gaussian of a=1.5 pixel (approx 0.7"). The circles 
are centered on the stacking position and have a radius of 2". The 
detection significance of the summed counts is 4.2cr for the “dusty” 
and < 2cr for the “old”. The bright spot in the right panel at about 4" 
from the stacked position, corresponds to an individually detected 
object. 

extreme than those inferred from the 2 Ms CDF-N spectral 
analysis (Vignali et al. 2002). Given the small area sampled 
and the spectroscopic incompleteness, the present ERO sam¬ 
ple is not suitable for a reliable estimate of the AGN fraction 
among the EROs population. However, the present low detec¬ 
tion rate of AGN among EROs (1/21) confirms that the bulk 
of the EROs population is not related to active phenomena 
(see Sect. 2.2). 

In order to constrain the average X-ray properties of individ¬ 
ually undetected EROs, we have applied the stacking tech¬ 
nique (Brandt et al. 2001; Nandra et al. 2002) separately for 
the two classes of objects. The samples consist of 12 “dusty” 
and 8 “passive” EROs, with average redshifts of z = 1.053 
and 2 = 1.145, respectively. We performed the detection in 
the three standard X-ray bands (0.5-8, 0.5-2 and 2-8 keV) 
and in the 1-5 keV band which roughly corresponds to the 
rest-frame hard X-ray band at the average redshift of the 
samples. We refer to Brusa et al. (2002b) for details on the 
X-ray data analysis. 

In Fig. 2 we show the summed images in the full (0.5-8 keV) 
band. The most important finding is that the two spectroscop¬ 
ically classes of EROs have different X-ray properties. 

3.1. An X-ray dichotomy 

The X-ray emission from “old” EROs remains undetected in 
all the considered energy bands, for a total effective exposure 
time of ~ 7.0 Ms. Assuming a thermal emission model with 
kT=l keV and solar metallicity, fully consistent with the av¬ 
erage X-ray spectrum of nearby elliptical galaxies (e.g., Pel¬ 
legrini 1999), the 3cr upper limit on the 0.5-2 keV luminosity 
is Lx < 10 41 erg s _1 . 

An excess of counts above the expected background level is 
clearly present in the “dusty” sample, for a total effective ex¬ 
posure time of 10 Ms. The signal is stronger in the full band 
(99.997% confidence level assuming a Poisson distribution. 


corresponding to about 4.2a) and it is still present, although 
at a slightly lower confidence level, in the soft and 1-5 keV 
band, while it is not statistically significant in the hard band. 
Assuming a T = 2 power-law spectrum plus Galactic ab¬ 
sorption (N h = 8 x 10 19 cm -2 ), the 3a upper limit on the av¬ 
erage band ratio (H/S<0.85, where H and S indicates the hard 
band and soft band counts, respectively), corresponds to an 
absorption column density lower than 10 22 cm -2 (see Fig. 4 
in Alexander et al. 2002). Given that the hard X-ray spectrum 
of nearby starburst galaxies is well-represented by power- 
law models with slopes in the range 1.8-E2.5 and low in¬ 
trinsic absorption (e.g., Ptak et al. 1999; Ranalli et al. 2002), 
we adopted T=2.1 to compute X-ray fluxes and luminosities. 
With this assumption, the stacked count rate in the 1-5 keV 
band corresponds to an average 2-10 keV rest-frame lumi¬ 
nosity of ~ 8 x 10 4 ° erg s _1 , at the mean redshift z = 1.053. 
The inferred, average X-ray luminosity is around one order 
of magnitude larger than that of normal spiral galaxies (Mat- 
sumoto et al. 1997), it is similar to that of the starburst galaxy 
M82 (e.g., Griffiths et al. 2000; Kaaret et al. 2001) and about 
3 times lower than that of NGC 3256, one of the most X-ray 
luminous starburst in the local Universe (e.g., Moran, Lehn- 
ert, & Helfand 1999). 

The starburst scenario is supported also by the optical data. 
Indeed, following Cimatti et al. (2002), we recomputed the 
average spectrum of the 12 EROs we used in the X-ray stack¬ 
ing analysis. The average optical spectrum shows a very red 
and smooth continuum, with a strong |OII] A3727 emission 
and lacks of the presence of clear AGN indicators such as 
[NeV]A3426, suggesting that the optical emission is domi¬ 
nated by star-forming systems. 

4. Hard X-ray properties of EROs 

We have collected from the literature a sample of hard X- 
ray selected EROs serendipitously detected in moderate-deep 
and deep Chandra and XMM -Newton observations (Horn- 
schemeier et al. 2001; Alexander et al. 2002; Mainieri et al. 
2002). From the observed X-ray band ratio, we have com¬ 
puted for each ERO the corresponding X-ray column den¬ 
sity for a source with a typical AGN power-law spectrum 
(T=1.8) at z — l, while for the sources in the Lockman Hole 
we adopted the best-fit value quoted by Mainieri et al (2002) 
and corrected to 2 = 1 for unidentified sources. The results 
are reported in Fig. 3. Almost all of the individually detected 
sources are consistent with intrinsic column densities in ex¬ 
cess of 10 22 cm -2 . The so far individually detected EROs 
have high X-ray luminosity (Lx > 1 0 13 erg s _1 , in a few case 
even larger than 10 44 erg s _1 ) and they actually are heav¬ 
ily obscured AGN, as inferred from X-ray spectral analysis 
(Cowie et al. 2001; Gandhi et al. 2002; Mainieri et al. 2002). 
Therefore, hard X-ray selected EROs (or at least a fraction of 
them) have properties similar to those of Quasars 2, the high- 
luminosity, high redshift type II AGNs predicted in X-Ray 
Background synthesis models (e.g. Comastri et al. 1995). 
Stacking analysis of those objects not individually detected 
in the X-rays suggests a softer X-ray spectrum than that of 
AGN-powered EROs (see also Alexander et al. 2002). More- 
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Fig. 3. X-ray absorption column density vs. hard X-ray flux for 
a sample of AGN-powered EROs, serendipitously detected in hard 
X-ray surveys. Different symbols refer to different samples as fol¬ 
lows: triangles = “Daddi” field; squares = CDFN (Homschemeier et 
al. 2001; Alexander et al. 2002); circles = Lockman Hole (Mainieri 
et al. 2002). The individually detected AGN in the K2 0 CDFS area 
is marked with a pentagon; the result from the stacking analysis of 
the “dusty” sample is marked with an asterisk. Dot-filled symbols 
refer to sources with measured redshift. 


over, the results from the CDFS K2 0 EROs sample suggest 
that different spectroscopic classes of EROs are characterized 
by different high-energies properties. Indeed, “dusty” EROs 
are moderately luminous hard X-ray emitters and their av¬ 
erage X-ray spectrum (N# < 10 22 cm -2 ) is consistent with 
that measured for nearby star-forming galaxies (e.g. Dahlem, 
Weaver & Heckman 1998). Early-type galaxies among EROs 
are not detected in the X-rays and their 3cr upper limit on the 
soft X-ray luminosity is consistent with the emission from 
nearby elliptical galaxies (Pellegrini 1999). 


5. Conclusions 

The most important results can be summarized as follows; 

• Hard X-ray observations turned out to be extremely use¬ 
ful to pick-up AGN-powered EROs in optically selected 
samples. The fraction of AGN among EROs is in the 
range 1 A 15%, depending on the limiting fluxes in the 
X-ray and near-infrared bands. 

• AGN-powered EROs are X-ray luminous and obscured 
objects. The intrinsic column densities inferred for a sam¬ 
ple of hard X-ray selected EROs are in excess of 10 22 
cm -2 . 

• “Dusty” EROs are relatively unobscured (N// < 10 22 
cm -2 ) hard X-ray sources and their X-ray properties are 
consistent with those of nearby star-forming galaxies. 

• “Old” EROs are not detected in the X-rays. 
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